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Abstract 

This paper aims to show that making use of Newton's view on equations of motion of a 

physical system and of the Maxwell stress tensor we come to a natural nonlinearization of 

Maxwell equations in vacuum making use only of nonrelativistic terms. The new equations 

include all Maxwell solutions plus new ones, among which one may find time-stable and 

»j . spatially finite ones with photon-like properties and behavior. 

GC 

■jh , 1 Introduction 

>y 

r£3 ' As it is well known the vacuum Maxwell equations (zero charge density: p = 0) do not admit 

spatially finite time-stable solutions of photon-like type. This is due to the fact that every com- 
ponent U(x,y,z,t) of the electric E and magnetic B fields necessarily satisfies the D'Alembert 
wave equation \DU = 0, and according to the Poisson's theorem for this equation, every spatially 
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finite initial condition U(x, y, z, 0) = ip(x, y, z); tsHcc, y, z, 0) = i[)(x, y, z), where ip and ip are fi- 



at 
nite functions, blows up radially and goes to infinity with the speed of light [1,2,3]. So, through 

every spatial point outside the initial condition pass fore-front and back-front, and after this 



the corresponding point forgets about what has happened. This rigorous mathematical result 
does not allow Maxwell vacuum equations to describe finite electromagnetic pulses propagating 
uniformly in vacuum as a whole along some spatial direction without dispersion. Moreover, no 
expectations for description of photon-like objects having besides translational also rotational 
"55 ■ component of propagation would seem to be reasonable and well-grounded. 

£>->! On the other hand the Poynting energy-momentum balance equation in vacuum 

where c is the velocity of light in vacuum, admits time-stable spatially finite solutions with 
3-dimensional soliton-like behavior, for example (no rotation component of propagation) 

E = [u(x, y, ct + sz), p (x, y, ct + ez), 0], B = [ep (x, y, ct + ez), —eu(x, y, ct + ez),0], e = ±1 , 

where u and p are arbitrary functions, so they can be chosen finite. This observation suggests 
to look deeper and more carefully into the structures and assumptions used for mathematical 
interpretation of the experimental electric-magnetic induction discoveries made in the 19th cen- 
tury. In other words, which relations and on what grounds should be defined as basic, so that 
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the further deduced equations and relations to give reasonable and physically acceptable results. 
Finding the right way to choosing adequate mathematical objects and corresponding equations 
seems specially important when we try to describe the intrinsic dynamical properties of such 
spatially finite time stable field objects. Therefore, it seems very important to have the right 
notion about concepts like physical object, intrinsic property, dynamical property, identifying 
characteristics, admissible changes, field equations, etc. Some preliminary considerations in this 
direction might be the following ones. 

From a definite point of view every physical system is characterized by two kinds of prop- 
erties. The first kind of properties we call identifying, they identify the system throughout its 
existence in time, so, the corresponding physical quantities/relations must show definite conser- 
vation/constancy properties (with respect to the identification procedure assumed). Without 
such experimentally established properties we could not talk about physical objects/systems 
at all. The second kind of properties (which may be called kinematical) characterize the 
time-evolution of the system, the corresponding quantities are time-dependent, and the cor- 
responding evolution is consistent with the conservative/constant character of the identifying 
properties/quantities. In this sense, the equations of motion of a physical system can be consid- 
ered as relations determining the admissible time-changes of these quantities. For example, the 
mass m of a classical particle is an identifying quantity, while the velocity v is a kinematical one. 
This view implies, of course, that the external world acts on the system under consideration also 
in an admissible way, i.e. an assumption is made that the system survives, the interaction with 
the outside world does not lead to its destruction. 

In theoretical physics we usually make use of quantities which are functions of the identifying 
and of the kinematical characteristics of the system and call them dynamical quantities. A 
well known example is the momentum p of a particle: p = mv. Of crucial importance for the 
description of admissible changes are the conservative dynamical quantities, i.e. those which may 
pass from one physical system to another with NO LOSS. For example energy and momentum 
are such quantities, moreover, they are universal in the sense that every physical object carries 
non-zero energy-momentum and, vice versa, every quantity of energy-momentum is carried by 
some physical object. So, if a definite quantity of energy-momentum passes from one object to 
another, this same quantity of energy-momentum can be expressed in terms of the characteristics 
of the two objects/systems, and the two expressions to be equalized. In this way we obtain 
consistent with the energy-momentum conservation law equations of motion, and this is the way 
used by Newton to write down his famous equations p = F, where F carries information about 
where the momentum change of the particle has gone, or has come from. This also clarifies 
the physical sense of the concept of force as a change of momentum, or as a change of energy- 
momentum in relativistic terms. Paying due respect to Newton we shall call some equations of 
motion of Newton type if on the two sides of " = " stay physical quantities of energy-momentum 
change, or energy-momentum density change in the case of continuous systems. Note that, 
written down for the vector field p, i.e. in terms of partial derivatives, the above Newton 
equation looks like V p p = rnF, where the left hand side means performing two steps: first, 
determining the "change quantity" Vp, second, projecting Vp on p, and the right hand side 
may be expressed as a function of the characteristics of both: the particle and the external 
physical envirenment. 

If there is no energy-momentum (or energy-momentum density) change, then putting the 
corresponding expression equal to zero, e.g. V p p = 0, we obtain the "free particle" or "free 
field" equations. In such a case we just declare that only those changes are admissible which are 
consistent with the (local and integral) energy-momentum conservation. 

We note that an initial extent of knowledge about the system we are going to describe 



mathematically is presupposed to be available, so that the assumptions made to be, more or 
less, well grounded. This knowledge is the base that generates corresponding insight and directs 
our attention to the appropriate mathematical structures. This is exclusively important when 
we deal with continuous, or field, physical objects/systems. 

In view of the above considerations, roughly speaking, in the free field case the steps to follow 
are: 

1. Specify and consider the mathematical model-object $ which is chosen to represent the 
integrity of the physical system considered; 

2. Define the change-object D(&); 

3. "Project" D(&) on $ by means of some (in most cases bilinear) map *P; 

4. The projection *p(D(<E>), <£) obtained we interpret physically as energy-momentum change; 

5. We put this projection equal to zero: *p(D($), <I>) = . 

The zero value of the projection *}3(Z)(3>), <E>) is interpreted in the sense that the identifying 
characteristics of <I> have not been disturbed, or, the change D(<&) is qualified as admissible. This 
consideration shows the importance of knowing how much and in what way(s) a given physical 
system is potentially able to lose, or gain energy-momentum (locally or globally), without losing 
its identity. 

It is always very important to take care of the physical sense of the quantities that we put 
on the two sides of the relation A = B. Mathematically, from set theory point of view, A and 
B denote the same element, which element may be expressed in different terms, e.g. the real 
number 2 can be expressed as 3 — 1 = 6/3 = 2(sin 2 x + cos 2 x) = ^(2x + const) and also in 
many other ways. From physical point of view, however, we must be sure that A and B denote 
the same thing qualitatively and quantitatively, i.e. the same physical quantity. This is specially 
important when the equation we want to write down constitutes some basic relation. And the 
point is not the physical dimension of the two sides to be the same: any two quantities by means 
of an appropriate constant can be made of the same physical dimension, but this is a formal 
step. The point is that the physical nature of the physical quantity on the two sides must be the 
same. 

For example, it is quite clear that on the two sides of the Newton's law p = F stays the well 
defined for any physical system quantity " change of momentum" since the momentum quantity 
is a universal one. For a counterexample, which physical quantity stays on the two sides of the 
Poisson equation AU = kp? On one hand, such a quantity is expressed through AU and, since 
gradU is usually interpreted as force, AU appears as a "change of force" characteristic of the 
field U since it is essentially defined by the second derivatives of U. On the other hand, the 
same quantity is expressed through kp and appears as a characteristic of the mass particles, so, 
do we know such a quantity? The same question can be raised for one of the Maxwell equations: 
rotB- ±E = ^j. 

c c •> 

In the case of classical particles momentum is always represented as the product mv and this 
is carried to fluid mechanics as p(x,y,z;t).v(x,y,z;t), where p is the mass density. A similar 
quantity is introduced in electrodynamics as electric current density j = p(x,y, z;t).v(x,y, z;t), 
where p is the electric charge density. The energy-momentum exchange between the field and the 
available charged particles is described by the force field F = /?E + -j x B. So, the corresponding 
Faraday-Maxwell force lines should be the integral lines of the vector field F. Clearly, in the 
charge-free case we get F = 0, so the concept of force-lines defined by F does not work. Hence, 
if we would like to use this concept appropriately in the charge-free case, we have to introduce 
it appropriately. The simplest way seems to consider the integral lines of E and B as force lines 
also in the charge free case, but we do not share this view: if p = then j = 0, the force-vector 



is zero and NO integral force lines exist. The two vectors E and B generate, of course, integral 
lines, but these integral lines are NOT force lines since E and B are NOT force fields, and such 
an interpretation of the integral lines of E and B would be misleading. In fact, there exist a 
sufficiently good force field defined by Maxwell in terms of the divergence of his stress tensor 
M !J , which definition works quite well also out of and away from any media built of, or 
containing, charged mass particles. So, in the frame of the theory at the end of 19th century 
if we ask the question: if there are NO charged particles and the time- dependent EM-field cannot 
transfer energy-momentum to them by means of the force field F, and the propagation of the 
free EM-field is available, so that energy-momentum internal exchanges should necessarily take 
place, how these processes and the entire propagational behaviour of the field could be understood 
and modeled?, the right answer in our view should be: turn to M %3 and consider carefully 
the divergence VjM 1 - 7 as corresponding force field generating corresponding force 
lines along which energy-momentum is locally transported. As will be seen further 
in the paper, such a look on the issue would necessarily lead Maxwell and his followers to the 
prediction that real, free, spatially finite and time-stable formations of electromagnetic field 
nature having translational-rotational dynamical structure should exist, a result that has been 
proved in studying the photoeffect phenomena about 30 years after Maxwell's death. 

We consider as a remarkable achievement of Maxwell the determination of the correct ex- 
pressions for the energy density of the electromagnetic field through the concept of stress [4] . His 
electromagnetic stress tensor M %3 still plays an essential role in modern electromagnetic theory 
as a part of the modern relativistic stress-energy-momentum tensor. However, by some reasons, 
Maxwell did not make further use of the computed by him divergence VjM 1 - 7 of the stress tensor 
(and called by him "force field" [4]) for writing down Newton type equations of motion for a free 
electromagnetic field through equalizing different expressions for the same momentum change. 
Probably, he had missed an appropriate interpretation of the vector cE x B (introduced by 
Poynting 5 years after Maxwell's death and called "electromagnetic energy flux" [5]). 

In this paper we consider one possible approach to come to natural free field equations of 
motion of Newton type that could be deduced making use of formally introduced vacuum analog 
of the Maxwell stress tensor and of the Poynting vector in the frame of the available theoretical 
notions and concepts at the end of 19th century. As a first step we are going to show that an 
analog of Maxwell's stress tensor participates in a (well known today) mathematical identity 
having nothing to do with any physics. 

2 A non-physical view on Maxwell stress tensor 

The mathematical identities have always attracted the attention of theorists, in particular, 
those identities which involve the derivatives of the objects of interest (differential identities). 
A well known such example is the Bianchi identity satisfied by any connection components: 
this identity is a second order system of (in general, nonlinear) partial differential equations. 
The gauge interpretation of classical Maxwell electrodynamics, as well as the Yang-Mills theory, 
substantially make use of this identity. Such identities are of particular importance when on the 
two sides of "=" stay correctly (i.e. in a coordinate free way) defined expressions. 

It is elementary to show that in the frame of classical vector analysis any two vector fields 
(V, W) and the corresponding mathematical analog of the Maxwell stress tensor M %d (V, W) are 
involved in a differential identity. Introducing the Maxwell stress tensor in such a formal way 
in the vacuum case will help us avoid all questions concerning the structure and properties of 
aether. 

We begin with the well known differential relation satisfied by every vector field V on the 



euclidean space M 3 related to the standard coordinates (x i = x, y, z), i = 1, 2, 3, denoting by V 2 
the euclidean square of V, by " x " - the vector product, and using the V-operator: 

- V(V 2 ) = V x rot V + (V.V)y = V x rot F + V V V. 

Clearly, on the two sides of this relation stay well defined quantities, i.e. quantities defined 
in a coordinate free way. The first term on the right hand side of this identity accounts for 
the rotational component of the change of V, and the second term accounts mainly for the 
translational component of the change of V. Making use of component notation we write down 
the last term on the right side as follows (summation over the repeated indices): 

(V v Vy = V i V i V j = ViiV^) - V'ViV* = Vi{V l V j ) - f'div V. 

Substituting into the first identity, and making some elementary transformations we obtain 

Vi (v*V j - h ij V 2 j = [(rot V) x V + VdivV} : 
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where b %3 = 1 for i = j, and S 13 = for i ^ j. If now W is another vector field it must satisfy 
the same above identity: 

Vi (w { W j - ]-5 ij W 2 J = [(rot W) x W + WdivW] j . 

Summing up these two identities we obtain the new identity 

/ V 2 + W 2 

ViM 11 = V. V l V 3 + W l W ] - b 13 ■ 



2 
- [(rot V) x V + VdwV + (rot W) x W + WdwW] J . (1) 

We emphasize once again the two moments: first, this identity (1) has nothing to do with 
any physics; second, on the two sides of (1) stay well defined coordinate free quantities. We 
note also the invariance of M %3 with respect to the transformations (V, W) — ► (—W,V) and 
(V,W)->(W,-V). 

The expression inside the round brackets on the left of (1), denoted by M 13 , looks formally 
the same as the introduced by Maxwell tensor from physical considerations concerned with the 
electromagnetic stress energy properties of continuous media in presence of external electromag- 
netic field. This allows to call formally any such tensor Maxwell stress tensor generated by 
the two vector fields (V, W) . The term "stress" in this general mathematical setting could be 
interpreted (or, justified) in the following way. Every vector field on M 3 generates corresponding 
flow by means of the trajectories started from some domain U C M 3 : at the moment t > 
the domain U is diffeomorphically transformed to a new domain Ut C M 3 . Having two vector 
fields on M 3 we obtain two consistent flows, so, the points of any domain U C M 3 are forced to 
accordingly move to new positions. 

Physically, we say that the corresponding physical medium that occupies the spatial region 
U and is parametrized by the points of the mathematical subregion U C M 3 , is subject to 
consistent and admissible physical "stresses" generated by physical interactions mathematically 
described by the couple of vector fields (V, W), and these physical stresses are quantitatively 
described by the corresponding physical interpretation of the tensor M l3 (V, W). 



We note that the stress tensor M 1J in (1) is subject to the divergence operator, and if we 
interpret the components of M %3 as physical stresses, then the left hand side of (1) acquires 
in general the physical interpretation of force density. Of course, in the static situation as it 
is given by relation (1), no energy-momentum propagation is possible, so at every point the 
forces mutually compensate: VjM 1 - 7 = 0. If propagation is allowed then the force field is NOT 
zero: VjM* J ^ 0, and we may identify the right hand side of (1) as a real time-change of 
appropriately defined momentum density P. So, assuming some expression for this momentum 
density P we are ready to write down corresponding field equation of motion of Newton type 
through equalizing the spatially directed force densities ViM 13 with the momentum density 
changes along the time coordinate, i.e. equalizing V 'jM !J with the ct-derivative of P, where 
c = const is the translational propagation velocity of the momentum density flow of the physical 
system (V, W). In order to find how to choose P we consider briefly the eigen properties of M %3 . 

3 Eigen properties of the Maxwell stress tensor 

We consider M lJ at some point of M 3 and assume that in general the vector fields E and B 
are linearly inependent, so E x B / 0. Let the coordinate system be chosen such that the 
coordinate plane (x,y) to coinside with the plane defined by E,B. In this coordinate system 
E = (Ei,E2,0) and B = (i3i,l? 2 ,0), so, identifying the contravariant and covariant indices 
through the Euclidean metric (so that M lJ = M % - = Mij), we obtain the following nonzero 
components of the stress tensor: 

M\ = (E 1 ) 2 + (B 1 ) 2 - -(E 2 + B 2 ); M\ = M 2 = E 1 E 2 + B x B 2 ; 

M 2 2 = (E 2 ) 2 + (B 2 ) 2 - -(E 2 + B 2 ); Ml = --(E 2 + B 2 ). 

Since M\ = — M|, the trace of M is Tr(M) = — ^(E 2 + B 2 ). The eigen value equation acquires 
the simple form [(M^) 2 — (A) 2 ] + (M 2 l ) 2 1 (M| — A) = 0. The corresponding eigen values are 

-Ly O - Ox * / f _ -1 x n / - -1 x r» . i 



X 1 = --(E 2 + B 2 ); A 2 , 3 = ±V( M i) 2 + (MlY = ±-y/{h) 2 + (/ 2 ) 2 , 

where I\ = B 2 — E 2 , / 2 = 2E.B. The corresponding to Ai eigen vector Z\ must satisfy the 
equation E(E.Zi) + B(B.Zi) = 0, hence, Z\ is proportional to E x B: Zi = fcExB. The other 
two eigen vectors Z\ 2 satisfy correspondingly the equations 



e(e.z 1i2 ) + b(b.z 1)2 ) = [ ± ^(i 1 y + {i 2 y + I(e 2 + b 2 ; 



Zl,2- 



Taking into account the relation 



,E 2 + B 2 \ 2 , l2 E 2 + B 2 , 
{h? + {h? = z - |E x B| 2 , so |E x B| > , 



1 r 
4 

we conclude that the coefficient before Z\^ on the right is always different from zero, therefore, 
the eigen vectors Z\^ lie in the plane defined by (E, B). 

The above consideration suggests that the intrinsically defined potential dynamical abilities 
of propagation of the field are: translational along (E x B) ; and rotational inside the plane 
defined by (E,B). 



We give some motivation now that the theory should choose A2 = A3 = 0. In fact, it turns 
out that if these two relations do not hold then the translational velocity of propagation is 
allowed to be less then the speed of light in vacuum c. Recall first the transformation laws of 
the electric and magnetic vectors under Lorentz transformation defined by the 3-velocity vector 



v and corresponding parameter j3 = v/c. If 7 denotes the factor 1/y 1 — /3 2 then we have 

E' = 7 E + -^P- v(E.v) + -v x B, 

V 1 c 

B' = 7B + -Z2 v (B.-v) - -V X E. 
V 1 c 

Assume first that I2 = 2E.B = 0, i.e. E and B are orthogonal, so, in general, in some 
coordinate system we shall have ExB/0. 

If Ii > 0, i.e. |E| < |B|. We shall show that the choice - = (|E|/|B|) < 1 is addmissible. 
Actually, the assumptions E' = and v.B = lead to 7 v.E + (1 — 7)(E.v) = 0, i.e. E.v = 0. 
Thus, E' = leads to c[E| = v|B||sin(v, E)|, and since v.E = then |sin(v, E)| = 1. It follows 
that the speed v = cfgr < c is allowed, so the vector product E' x B' = 0. 

If I\ < 0, i.e. |E| > |B|, then the choice B' = and v.E = analogically leads to the 

l"R I 

conclusion that the speed v = cWf < c is allowed and E' x B' = 0. 

Assume now that I2 = 2E.B 7^ 0. We are looking for a reference frame K' such that 
E' x B' = 0, while in the reference frame K we have E x B ^ 0. We choose the relative velocity 
v such that v.E = v.B = 0. Under these conditions the equation E' x B' = reduces to 

E xB + -(E 2 + B 2 ) = 0, so, - = |Ex B|/(E 2 + B 2 ). 
c c 

Now, from the above mentioned inequality E 2 + B 2 — 2|E x B| > it follows that - < 1. 

These considerations show that under nonzero I\ and I2 the translational velocity of propa- 
gation of the field, and of the energy density of course, will NOT be equal to c. Hence, the only 
realistic choice left is I\ = I2 = 0, which is equivalent to E 2 + B 2 = 2|E x B|. Hence, assuming 
|Tr(M)| to be the energy density of the field, the name "electromagnetic energy flux" for the 
quantity cE x B seems well justified without turning to any field equations. 

These considerations show also that if ii = |B| 2 — |E| 2 = then the electric and magnetic 
components of the field carry always the same energy density, so, a local mutual energy exchange 
is not forbidden in general, but, if it takes place, it must be simultanious and in equal quantities. 
Hence, under zero invariants I± = and I2 = 2E.B = internal energy exchange is allowed but 
this exchange occurs without available interaction energy, and the field equations must take care 
of such a possibility. 

4 Nonlinear equations for the electromagnetic field 

We procede to write down dynamical equations of the field through specializing how the internal 
local momentum exchange is realized. 

We note first that the assumption that the energy density of the field coincides with \Tr (M) \ = 
i [E 2 + B 2 ] presupposes that there is NO interaction energy between the electric and magnetic 
components of the field. It is an usual practice to introduce interaction energy through some 
function of the scalar product of the two fields: Ei n t = /(E.B), e.g., proportional to E.B. Hence, 
in the free field case it seems natural to assume E.B = 0, i.e. I2 = 0. We just note that, this 
does not mean that there is no energy exchange between the electric and magnetic components. 



Now, what about a possible momentum exchange between the electric and magnetic compo- 
nents of the field? Note that the assumption that the field momentum is given by -E x B, i.e. 
it is a bilinear function of the electric and magnetic components, clearly suggests that there is 
available nonzero interaction momentum. The point now is to get some clarification how this 
local momentum exchange takes place and to find appropriate mathematical representatives of 
the corresponding partners realizing such an exchange since neither E nor B are able to carry 
momentum separately. 

We procede to answer this question through writing down corresponding field equations. 

Replacing (V, W) in (1) with (E, B) we obtain 

ViAT J = Vi ( E*E^' + B^ - S ij 



2 
- [(rot E) x E + Ediv E + (rot B) x B + Bdiv B] 3 . (2) 

As we mentioned, in the static case, i.e. when the vector fields (E, B) do not depend on the 
time coordinate £ = ct, NO propagation of field momentum density P should take place, so, 
at every point, where (E, B) ^ 0, the stress generated forces must mutually compensate, i.e. 
the divergence VjM lJ should be equal to zero: VjM 1 - 7 = 0. In this static case Maxwell vacuum 
equations 

rotE + — = 0, rotB- — = 0, divE = 0, divB = (*) 

give: rotE = rotB = 0; divE = divB = 0, so, all static solutions to Maxwell equations determine 
a sufficient, but NOT necessary, condition that brings to zero the right hand side of (2) through 
forcing each of the four vectors there to get zero values. 

In the non-static case, i.e. when 4r- ^ 0; ^- 7^ 0, time change and propagation of field 
momentum density should take place, so, a full mutual compensation of the generated by the 
Maxwell stresses at every spatial point local forces may NOT be possible, which means VjM^ ^ 
in general. These local forces generate time-dependent momentum propagation P(E, B) at the 
spatial points. Therefore, if we want to describe this physical process of field momentum density 
time change and spatial propagation we have to introduce explicitly the dependence P(E,B). 
If we follow the classical (nonrelativistic) way of consideration and denote by J the vector field 
with components & = V 'jM iJ , we can write down the force flow across some finite 2-surface S 
in the usual (and widely spread in almost all textbooks) way as f s $.ds (from modern point 
of view we should write ig(dx A dy A dz) instead of 5-ds under the integral, where i$ denotes 
the inner product between the vector field 5 and the volume form dx A dy A dz, i.e. to make 
use of the Poincare isomorphism between vector fields and 2-forms on R 3 ). This flow generates 
changes of the momentum density flow across S which should be equal to j| J s P(E, B).ds. We 
obtain 

j t /p(E,B).ds = J Z.ds. 

The explicit expression for P(E, B), paying due respect to J.Poynting [5], and to J.J.Thomson, 
H. Poincare, M. Abraham [6], and in view of the more than a century development of the elec- 
tricity theory and practice, has to be introduced by the following 

Assumption: The field momentum density is given by P := -E x B . 



According to the Assumption and the above interpretation of the relation VjM JJ 7^ 0, and 
in view of the arbitrariness of the 2-surface S we come to the vector differential equation 

|(ExB)=5, e = ct, (**) 

which according to relation (2) is equivalent to 

rotE + — I x E + EdivE+ (rotB- — | x B + BdivB = 0. (3) 

This last equation (3) we write down in the following equivalent way: 



f)~pt 
rotE + — ) x E + BdivB 



dE , 
rotB - — ) x B + EdivE 



(4) 



The above relation (**) and the corresponding differential relation (3)/(4) we consider as math- 
ematical adequate in momentum-change terms of the electric-magnetic and magnetic-electric 
induction phenomena in the charge free case. We recall that these induction phenomena are 
described in what we call "Faraday-Maxwell theory" by the following well known integral and 
differential equations 



l/ s B - 



rotE.ds 



<9B 



— — = — rotE, (the Faraday induction law), 



rotB, (the Maxwell displacement current law). 



d [^ , f „ , dE 

— / E.ds = / rotB.ds ->• — - 

<% Js Js d£ 

We stress once again that these last Faraday-Maxwell relations have NO direct energy-momentum 
change-propagation (i.e. force flow) nature, so they could not be experimentally verified in a 
direct way. Our feeling is that, in fact, they are stronger than needed. So, on the corresponding 
solutions of these equations we'll be able to write down formally adequate energy-momentum 
change expressions, but the consistency of these expressions with the experiment will crucially 
depend on the nature of these solutions. As we already mentioned, the nature of the free solu- 
tions (with no boundary conditions) to Maxwell vacuum equations with spatially finite initial 
conditions requires strong time-instability (the Poisson theorem for the D'Alembert wave equa- 
tion). And time-stability of time-dependent vacuum solutions usually requires spatial infinity 
(plane waves), which is physically senseless. Making calculations with spatially finite parts of 
these spatially infinite solutions may be practically acceptable, but from theoretical viewpoint 
assuming these equations for basic ones seems not acceptable since the relation "time stable 
physical object - exact free solution" is strongly violated. 

Before to go further we write down the right hand side bracket expression of (4) in the 
following two equivalent ways: 



rotB + 



0(-E) 



xB + (-E)div(-E) 



rot 1 






-B) + EdivE 



(5) 



These last two expressions (5) can be considered as obtained from the left hand side of (4) under 
the substitutions (E,B) — > (B, — E) and (E,B) — > (— B,E) respectively. Hence, the field (E,B) 
has always as a partner one of the fields (— B, E), or (B, — E). 

This observation suggests the following view: an adequate mathematical representation of a 
time dependent free electromagnetic field requires a collection of two fields : (E, B); (— B, E) , or 



(E, B); (B, — E) . We could also say that a real free field consists of two interacting subsystems 

(partners) described by two component-fields, and each component-field has electric and magnetic 
components, and is determined by the other through rotation-like transformation, in particular, 
if E.B = then this rotation is to ±7r/2. This view and relation (4) suggest, in turn, that 
the intrinsic dynamics of the real time- dependent electromagnetic fields could be considered as 
accompanied by a local energy-momentum exchange between the corresponding two component- 
fields. Such a view suggests also that each of the two component-fields of a real free field may 
keep locally its energy-momentum if the inter-exchange is simultaneous and in equal quantities. 
We are going now to interpret the equation (4) in accordance with the view on equations 
of motion as stated in the Introduction. Our object of interest <3?, representing the integrity of 



a real electromagnetic field, is the couple 



(the other case 



(E,B);(B,-E) 



(E,B);(-B,E 

is considered analogically). In view of the above considerations our equations should be able 
to express first, internal for each component-field energy-momentum redistribution, second, 
possible and admissible energy-momentum exchange between the two component-fields. Hence, 
we have to define the corresponding change-objects D(E, B) and D{— B, E) for each component- 
field and their "projections" on the two component-fields. 

The change object D(E, B) for the first component-field (E, B) we, naturally, define as 

D(E,B) := (rotE + ^; divBJ . 



The corresponding "projection" of D(E, B) on (E,B) 



dB 



(rotE+— ;divBJ;(E,B) 



(BB \ 

rotE + — ;divBj;(E,B) 



dB\ 
rot E + — - x E + Bdiv B. 



<P[D(E,B);(E,B)]=<P 
is suggested by the left hand side of (4) and we define it by : 

For the second component-field (— B,E), following the same procedure we obtain: 
<PpD(-B,E);(-B,E)]=<P I rot(-B) + — ; divEJ ; (-B,E) = 

rot(-B) + — 1 x (-B)+EdivE= (rotB- — J x B + EdivE. 

Hence, relation (4) looks like 

<£ [£>(E, B); (E, B)] + % [D(-B, E); (-B, E)] = 0. 

The accepted two-component view on a real time dependent electromagnetic field allows 
in principle admissible energy-momentum exchange with the outside world through any of the 
two component-fields. Hence, the above calculations suggest to interpret the two sides of (4) 
as momentum quantities that each component-field (E,B), or (— B,E), is potentially able to 
give to some other physical object and these quantities are expressed in terms of E, B and their 
derivatives only. In the case of free field, since no momentum is lost by the field, there are 
two possibilities: each component-field to keep the energy-momentum it carries, or one of the 
component-fields to change its energy-momentum at the expense of the other. If we denote 
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by An an d by A22 the allowed energy-momentum changes of the two component fields, by 
A12 the energy-momentum that the first component-field receives from the second component- 
field, and by A21 the energy-momentum that the second component-field receives from the first 
component-field, then according to the energy-momentum local conservation law we may write 
the following equations: 

An = Ai2 + A 2 i; A 22 = - (A 2 i + A12) , 

which is in accordance with the equation (4): An + A22 = 0. 

We determine now how the mutual momentum exchange between the two component-fields 
P(E,B) «-* P(-B,E)) or > P(E,B) <-* P(B,-E) i s performed, i.e. the explicit expressions for A12 and 
A21. The formal expressions are easy to obtain. In fact, in the case P(e,b) — * P(_be)i i- e - the 
quantity A21, we have to "project" the change object for the second component-field 

rot(-B) + — ; divEj 

on the first component-field (E,B). We obtain: 

rot(-B) + — I x E + BdivE = - I rotB- — 1 x E + BdivE . (6) 

In the reverse case P(_b,e) — * P(e,b)> i- e - the case A 12, we have to project the change-object 
for the first component-field given by 

£>(E,B) := f rotE + ^; divB j 

on the second component-field (— B,E). We obtain 

(/9B \ / <9B \ 

rotE + — J x (-B)+EdivB = - (rotE+ — J xB + EdivB. (7) 

So, the internal local momentum balance is governed by the equations 



<9B \ / <9B \ / r?E N 

rotE+— ) xE + BdivB = - f rotE + — ) xB + EdivB- (rotB- — ) xE + BdivE, 

(8) 



di J \ di J \ d£ 



rotB- — ) xB+EdivE= f rot B - — ) xE-BdivE+ I rot E + — j xB-EdivB. (9) 

These two vector equations (8)-(9) we consider as natural Newton type field equations. Ac- 
cording to them the intrinsic dynamics of a free electromagnetic field is described by two couples 
of vector fields, [(E, B); (— B, E)], or [(E, B); (B, — E)], and this intrinsic dynamics could be in- 
terpreted as a direct momentum exchange between two well defined subsystems mathematically 
described by these two component-fields. 

A further natural specilization of the above two vector equations (8)- (9) could be made if 
we assume that this internal momentum exchange realizes a dynamical equilibrium between the 
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two component fields, i.e. each component-field conserves its momentum : An = A22 = 0. In 
such a situation each component-field loses as much as it gains, so, equations (8)-(9) reduce to 

An = (rotE + — - ) x E + BdivB = 0, (10) 



rotB- — I x B + EdivE = 0, (11) 

(9B \ / r)E \ 

rotE + — ) x B-EdivB+ I rotB- — - ) x E-BdivE = 0. (12) 

di J V d U 

Equation (12) fixes, namely, that the exchange of momentum density between the two component- 
fields is simultanious and in equal quantities, i.e. a permanent dynamical equilibrium 
between the two component-fields holds: P(e,b) *-* P(-b,e)> or > P(e,b) +-*■ P(b,-e)- 

Note that, if equations (10) and (11) may be considered as field-equivalents to the zero force 
field (eqn. (11)) and its dual (eqn. (10)), this double-field viewpoint and the corresponding 
mutual energy-momentum exchange described by equation (12) are essentially new moments. 

Equations (10)-(12) also suggest that the corresponding fields are able to exchange energy- 
momentum with other physical systems in three ways. If such an exchange has been accom- 
plished, then the exchanged energy-momentum quantities can be given in terms of the charac- 
teristics of the other physical system (or in terms of the characteristics of the both systems, e.g. 
the Lorentz force pE + A x B), and to be correspondingly equalized to the left hand sides of 
equations (10)-(12) in accordance with the local energy-momentum conservation law. 

Finally, we give the 4-dimensional relativistic picture (details see in [7]). If the Minkowski 
pseudometric r\ has signature (—,—,—,+) and F^ = E*, F12 = B 3 , F13 = — B 2 , F23 = B 1 , 
and {*F) a p = — \e a p^ u F f11 ' , d is the exterior derivative, 8 = *d* is the coderivative, then the 
Maxwell stress tensor and its divergence are extended respectively to 



rpV 



1 



F^F va + (*F)^(*FY° and V V T» = F^SFf + (**V((J * F ) 



These expressions clearly and respectfully show the two-component (F, *i ? )-structure of the field. 
Equations (8)- (9) are extended correspondingly to 



* 



*V((J * Ff = - [F^(5 * FY + (*F)^(SF) 1 

F^SFf = F^ u (5 * FY + (*F)^(5FY- 
Equations (10)-(12) are extended correspondingly to 

F a ^(dF) a p li = (*F) lu/ (S*FY = 0, (*FY fi (d*F) a ^ = F^(8FY = 0, a < (3; 

(*F)°^(dF) a ^ + F^id * F) a ^ = (5 * FYF Ufl + (SF) v (*F) vli = 0, a < (3. 

5 Some Properties of the nonlinear solutions 

Clearly, all solutions to Maxwell pure field equations (*) are solutions to our nonlinear equations 
(8)-(9) and (10)-(12), we shall call these solutions linear, and will not be interested of them. In 
this section we shall concentrate on those solutions of (10)-(12) which satisfy the conditions 

<9B <9E 

rotE + — ^0, rotB-— - ^ 0, divE / 0, divB^O. 
ot, o£ 

12 



These solutions we call further nonlinear. We note some of the properties they have. 

1. E.B = 0; 

2. (rotE + H| J .B = 0; ( rotB — |j| J .E = 0, From these two relations the classical 
Poynting energy-momentum balance equation follows. 

The above two properties are obvious from equations (10) and (11). 

3. If (E, B) defines a solution then (E', B') = (oE - bB; 6E + aB), where a, b G R, defines 
also a solution. This property is immediately verified through substitution. 

4. E 2 = B 2 . To prove this, we first multiply equation (10) on the left by E and equation 
(11) by B (scalar products). Then we make use of the above properties 1 and 2 and of the vector 
algebra relation X.(Y x Z) = Z.(X xY). 

Properties (1) and (4) say that all nonlinear solutions to (10)-(12) are null fields, i.e. the 
two well known relativistic invariants I\ = B 2 — E 2 and I2 = 2E.B of the field are zero. 

5. B. [rot B - ||) - E. [rot E + f|) = B.rotB - E.rotE = 0. 

To prove this property we first multiply (vector product) (10) from the right by E, recall property 
1, then multiply (scalar product) from the left by E, recall again E.B = 0, then multiply from 
the right (scalar product) by B and recall property 4. 

Property (5) suggests the following consideration. If V is an arbitrary vector field on R 3 then 
the quantity V.rotV is known as local helicity and its integral over the whole region occupied by 
V is known as integral helicity, or just as helicity of V. Hence, property 5 says that the electric 
and magnetic components of a nonlinear solution generate the same helicities. If we consider 
(through the euclidean metric) E as 1-form on R 3 and denote by d the exterior derivative, 
then E A dE = E.rotE dx A dy A dz, so, the zero helicity says that the 1-form E defines a 
completely integrable Pfaff system. The nonzero helicity says that the 1-form E defines non- 
integrable Id Pfaff system, so the nonzero helicity defines corresponding curvature. Therefore 
the equality between the E-helicity and the B-helicity suggests to consider the corresponding 
integral helicities J R3 E A dE = J K3 B A dB (when they take finite nonzero values) as a measure 
of the spin properties of the solution. 



6. Example of nonlinear solution: 
E = 



z z 

x, y, ct + ez)cos(—K- — \- const), </>(x, y, ct + sz)sm(— k- — h const), 



B 



z z 

E(f>(x, y, ct + ez) sin(— k- — h const), —ecf)(x, y, ct + ez)cos(— k- — h const), 



where 4>(x, y, ct + ez) is an arbitrary positive function, l is an arbitrary positive constant with 
physical dimension of length, and e and n take values ±1 independently. The form of this solution 
shows that the initial condition is determined entirely by the choice of (j), and it suggests also to 
choose the initial condition 4>t=o(x, y, ez) in the following way. Let for z = the initial condition 
4>t=o(x, y, 0) be located on a disk D = D(x, y; a, b; r ) of small radius r , the center of the disk to 
have coordinates (a, b), and the value of 4>t=o(x, y, 0) to be proportional to the distance R(x, y, 0) 
between the origin of the coordinate system and the point (x, y, 0), so, R(x,y, 0) = \J x 2 + y 2 , 
and D is defined by D = {(x, y)|y / (x — a) 2 + (y — b) 2 < r }. Also, let 6d be the smoothed out 
characteristic function of the disk D, i.e. 9d = 1 everywhere on D except a very thin hoop-like 
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zone Bd C D close to the boundary of D where 0d rapidly goes from 1 to zero (in a smooth 
way), and Op = outside D. Let also the dependence on z to be given by be the corresponding 
characteristic function 9(z; 2tt1 ) of an interval (z, z + 2tt1 ) of length 2irl on the z-axis. If 7 is 
the proportionality coefficient we obtain 

4>(x, y, z, ct + ez) = j.R(x, y, O).0£>.6(ct + ez; 2irl ). 

We see that because of the available sine and cosine factors in the solution, the initial condition 
for the solution will occupy a helical cylinder of height 2irl , having internal radius of r and 
wraped up around the z-axis. Also, its center will always be R(a, b, 0)-distant from the z-axis. 
Hence, the solution will propagate translationally - along the coordinate z with the velocity 
c, and rotationally - inside the corresponding infinitely long helical cylinder because of the in- 
dependence of the available periodical multiples. The curvature K and the torsion T of the 
screwline through the point (x, y, 0) € D will be 

K _ R(x,y,0) m _ k1 



i? 2 0r,y,O) + /2' R?(x,y,0)+If 

The rotational frequency v will be v = c/2tt1 , so we can introduce the period T = 1/v and 
elementary action h = E.T, where E is the (obviously finite) integral energy of the solution 
defined as 3d-integral of the energy density (E 2 + B 2 )/2 = <j) 2 (see the figures on p. 62 in hep- 
tfi/0403244). 

This example presents also a completely integrable differential (2-dimensional) system, i.e. 
there exist two functions / and g such that the Lie bracket [E, B] can be represented in the 

27r/ 2 27r/ 2 

form /E + gB. The appropriately normed local helicities — -^E A dE = -^B A dB generate 
the integral helicity h = E.T, i.e. the elementary action, where E is the integral energy of the 
solution, T = 2irl /c, and, clearly, h = const. If we interpret h as the Planck's constant then the 
relation h = E.T is equivalent to the Planck's relation E = hv, and h appears as the (integral) 
spin of the solution. 

6 Discussion and Conclusion 

The main idea of the paper is that carrying out the Newton way for writing down dynamical 
equations for particles in mechanics to writing down dynamical equations for continuous field 
systems should naturally result to nonlinear partial differential equations even in non-relativistic 
theories. Moreover, clarifying the sense of the information included in these dynamical equations 
according to the Newton approach, we come to the conclusion formulated in the Introduction, 
namely, we have to mathematically describe those changes of the object considered which are 
qualified as admissible and consistent with the system's identification and with the local energy- 
momentum balance relations. In the case of "free" systems these relations represent the local 
energy-momentum conservation properties of the system. The energy- momentum characteristics 
are chosen because of their two important properties: they are physically universal and conser- 
vative. This means that every physical object carries nonzero energy-momentum and, vice versa, 
every quantity of energy-momentum is carried by some physical object. Also, if a physical ob- 
ject loses/gains some quantity of energy-momentum then some other physical object necessarily 
gains/loses the same quantity of energy-momentum. If this viewpoint is assumed, then the prob- 
lem of finding appropriate dynamical equations for an object reduces mainly to: first, getting 
knowledge of the potential abilities of the object considered to lose and gain energy-momentum; 
second, to create adequate mathematical quantities describing locally these abilities. 
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The electromagnetic field, considered as a continuous physical object of special kind, gives a 
good example in this direction since, thanks to Maxwell's fundamental and summarizing works, 
all the information needed is available. The notices of Poynting [5], and Thomson, Poincarc 
and Abraham [6], showing the importance of the (deduced from Maxwell equations) vector 
-E x B from local energy-momentum propagation point of view, has completed the resource 
of adequate and appropriate mathematical objects since it appears as natural complement of 
Maxwell stress tensor, and allows to write down dynamical field equations having direct local 
energy-momentum balance sense. However, looking back in time, we see that this viewpoint for 
writing down field equations has been neglected, theorists have paid more respect and attention 
to the "linear part" of Maxwell theory, enjoying, for example, the exact but not realistic, and 
even physically senseless in many respects, plane wave solutions in the pure field case. 

Therefore, not so long after the appearance of Maxwell equations the photoeffect exper- 
iments showed the nonadequateness of the linear part of Maxwell theory as a mathematical 
model of electromagnetic fields producing realistic model-solutions of free time-dependent fields. 
Although the almost a century long time development of standard quantum and relativistic 
quantum theories that followed, a reasonable model-solutions describing individual photons, 
considered as basic, spatially finite and time-stable objects, these theories have not presented so 
far. Nobody doubts nowadays that photons really exist, and this very fact suggests to try first 
classical field approach in finding equations admitting 3d-finite and time stable solutions with 
appropriate properties. 

The historical perspective suggests to follow the 4-potential approach, but modern knowledge 
and experience, and even the Maxwell stress tensor achievements, suggest some different views. 
In fact, we have all reasons to consider the microobjects as real as all other physical objects, 
so, no point-like charges and infinite field model-solutions should be considered as adequate. 
Since the 4-potential approach in Maxwell theory does not allow spatially finite and time stable 
pure field solutions with photon-like structure and behavior its interpretation as a basic concept 
does not seem to be appreciable. Also, the 4-potential approach excludes many solutions of the 
charge free Maxwell equations. For example, in relativistic terms the Coulomb field is given by 
the 2-form F = \dr A d£, dF = 0, its Minkowski-dual is *F = q sin 9 dO A dtp, d*F = 0, where 
F has a global 4-potential, but *F has NO global 4-potential. Now, the 2-parameter family of 
2-forms (J, *#) = (aF — b * F; bF + a * F),a,b S K, gives an infinite number of solutions to 
Maxwell equations d^ = 0, d * J = admitting NO global 4-potential. This suggests the view 
that the 4-potential can be used as a working tool (wherever it causes no controversies) but not 
as a basic concept. 

In conclusion, paying due respect to the Newton view on dynamical equations and to the local 
energy-momentum conservation law we based our approach on the Maxwell stress tensor and on 
the Poynting vector as natural quantities carrying the physically meaningful energy-momentum 
characteristics of the electromagnetic field. The natural description in these terms is based on 
two component-fields: [(E,B)], [(— B,E)], or [(E,B)], [(B, — E)], and mutual energy-momentum 
exchange between the two component-fields (in equal quantities) is considered as possible and 
explicitly accounted. The equations obtained include all solutions to the charge free Maxwell 
equations. A basic part of the new (nonlinear) solutions have zero invariants. Among these 
zero-invariant new solutions there are time-stable and spatially finite ones having photon-like 
properties and behavior. An analog of the Planck relation E = hv holds for these solutions, 
where the constant h appears as an integral helicity of such a solution. 

This study was partially supported by Contract (j> 15 15 with the Bulgarian National Fund 
"Science Research". 
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